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Abstract 

VT invariance of complex potentials Vix) = [V{—x)]* combines their real symmetry 
with imaginary antisymmetry. We describe a new exactly solvable model of this 
type. Its spectrum proves real, discrete and "three- fold" , e = = (2n + ajY, 

n = 0, 1, . . . , with ttj > and j = 1, 2, 3. 
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Recently, Daniel Bessis and Bender and Boettcher [0 proposed a certain 
modification of the interpretation of the one- dimensional Schrodinger-like differential 
equations. In a way inspired by the methodical relevance of the spatially symmetric 
harmonic oscillator in field theory they tried to weaken the standard assumption 
of hermiticity of the Hamiltonian H = H~^, i.e., in the strict mathematical language, 
the requirement of the essentially self-adjoint character of this operator. 

The new idea has found its first physical applications and immediate tests in field 
theory (cf., e.g., refs. [Q for explicit illustrations). Its use in the time- independent 
quantum mechanics has simultaneously been proposed 0. As long as the time- 
reversal operator T performs the mere Hermitian conjugation in the latter case, 
a generalization of bound states with definite parity {Vip = ±ip G L2) has been 
found in normalizable states with a parity-plus-time-reversal combined symmetry 
(VTi: = ±ije L2). 

In the PT— symmetric quantum mechanics a stability hypothesis Im E = has 
been tested semi-classically ||, numerically [jl], 0, analytically [0, |^ and perturba- 



tively 0, 0. One may even return to the related older literature, say, on the cubic 
anharmonic oscillator V{x) = uj"^ + igx^, all the resonant energies of which remain 



real and safely bounded below |TOl O. The similar complete suppression of decay 



has now been observed and/or proved in several other analytic models [jT2|, |T3|. 

Serious problems arise for the non-analytic PT— symmetric interactions [Q. An 
indirect clarification of the difficulty is being sought, first of all, in the various analytic 
PT— analogs of the usual square well ||, ^, |n| . In the present note, another exactly 
solvable example of such a type will be proposed and analyzed in detail, therefore. 

Let us start by characterizing solvable potentials by their so called shape invari- 
ance. The review lists nine of these models. Once we restrict our attention to 
the mere confluent-hypergeometric-type equations on full line, we are left with the 
spatially symmetric harmonic oscillator V(^*^)(r) = cj^ (r — h/uj)^ — u and with its 
Morse-oscillator partner 

y{Morse)^^~^ = + 5^ exp(-2/5r) - 2B{A - (3/2) exp(-/3r). (1) 
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The PT— symmetrization of the former model V^^^^ has been recalled in the pio- 
neering paper 0]. The latter case (P is to be considered here. 

The real and Hermitian version of the model is of a non-ceasine; interest 

in the current literature ||T6|. Its consequent PT— symmetrization necessitates a 
few preliminary considerations. In the first step, let us make a detour and return 
once more to the real and Hermitian V^^^^ and to its generalized three-dimensional 
S chrodinger-equat ion present at ion 



1 

i>{r) = E^{r), £=-- + «, a >0. (2) 



We may remind the reader that one can safely return back by choosing the trivial 
£ = — 1 and/or £ = 0. Still, beyond this special "full-line" case with parity (— l)^"*"^, 
the singularity at r = remains prohibitively strong for any a > 0. 



Fortunately, according to our recent letter |T^, an unexpected remedy may be 
found in the PT— symmetrization. It renders even the singular equation (0) tractable 
on the full line. The core of its PT— regularization lies in an analytic continuation in 
the coordinate r. This enables us to work in the complex plane with an appropriately 
chosen cut. The cut starts at the singularity in r = and may be directed, con- 
veniently, upwards. In such a setting one may integrate equation (0) along the real 
line which is only slightly shifted downwards. In the notation with r = s — ic, c > 
and with a new real quasi-coordinate s G {~oo, oo) the physical boundary conditions 
ip{—ic ± oo) = remain the same as before. 

A non-numerical key to the complexified one-dimensional problem lies in a 
surviving exhaustive solvability of its differential equation in terms of the confluent 
hypergeometric functions, 

e--'/2^(x) = Cir-"+^/2 ^F,[{2-2a~E/uj)/4,l-a;ujr'') + 

+ C2 r"+'/' iFi ((2 + 2a- E/oj)/A, 1 + «; r^) . (3) 

For the large coordinates r which lie within the wedges |Imr| < |Rer|, both these 
hypergeometric series grow as exp(c<jr^/2) or, alternatively, degenerate to a poly- 
nomial ||T^. Discarding the former possibility as manifestly violating our physical 
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boundary conditions we arrive at a slightly unusual formula for the binding energies 

E = Egn = uj{4:n + 2-2qa), q = ±1, n = 0,l,2,... 

and get also wave functions containing Laguerre polynomials, 
^(r) = const. r-''"+^/2g-c.rV2 ^(-ga) 



The spectrum remains real, discrete and bounded below and it reproduces the cur- 
rent, equidistant harmonic energies at a = 1/2 ||17|| . 

A change of the coordinates r = const x exp (const x) mediates a transition to the 
Morse forces v^^^"^^^\ We have to convert the asymptotically normalizable harmonic 
wavefunctions ipn,±{r) ~ exp(— r^/2), Rer — > ±00 to their Morse asymptotically 
normalizable counterparts ip{x). Using the explicit rule 



-I e 



i>{r) = y/r(p{x) (4) 



in the harmonic Schrodinger equation (0) we arrive at the PT— symmetrized Morse 
bound-state problem 

^2 



— r — 00 exp(4ix) — D exp{2ix) 
dx'' 



ip{x) = Eifix). (5) 



Our change of variables maps the straight integration contour r = s — ic with the real 
s G (—00, 00) and with the positive constant c > onto the deformed, down-bent 
curve 

C = {x = V — iu II V E {—7i/2, 7r/2), u = u{v) = ln(c/ cost;)} . (6) 

The circles centered at the origin are mapped upon segments of horizontal lines. The 
upward-running cut from r = to r = ioo is represented, due to the periodicity of 
the exponential, by an infinite family of the vertical lines. The whole plane of r is 
mapped on the single vertical strip in x, and the multi-sheeted Riemann surface in 
r is projected upon the whole complex x— plane. 

Our manifestly VT symmetric equation (|^) contains the Morse potential (|l]) 
where f3 = —2i, B = iuj and A = i{l — E/2uj). The old energy E is to be re- 
interpreted as a new coupling D = E. The old centrifugal- like parameter a = £ + 1/2 
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re-appears as a new momentum and specifies the Morse energy e = cP' . The original 
harmonic oscillator recipe a ^ E^'y^ replaced by the new equivalent rule D —* a^. 
This is because the new integration contour C preserves the asymptotic growth or 
decrease of the general solutions (|^). 

At first sight, the latter conclusion may seem slightly confusing. Its clarification is 
easy. In equation (j^) the real parameters v and u are such that the values of v remain 
bounded. The variation of v parameterizes not only s = c tanv but also both the two 
semi-infinite branches of u{v) G (Inc, oo). By construction, the harmonic-oscillator 
states with the even and odd quasi-parity g = ±1 and energies E = = 
a;(4n + 2 =F 2a), n = 0,1, . . . are mapped upon the new Morse bound states with the 
energies 

e = e^^\D) = {2m+lTD/2ujy, m = 0, 1, . . . . (7) 

Mutatis mutandis, also the wave functions are easily obtained. One only has to pay 
attention to the unavoided [|1^ crossings at the integers D/2u. 

Phenomeno logically, the new spectrum is fairly rich. For D/Au = M+a — 1/2 > 
with a small a G (0, 1) and an appropriate integer M > the set of the quasi-even 
energies splits in the two subsequences. The first one is finite, 

^mU-iP)/4 =ik + af, k = 0,l,...,M-l (8) 

and decreases with the number of nodes in fix). The infinite rest is increasing, 

e^}^,iD)/A=ik + l-af, k = 0,l,... . (9) 

In contrast, all the quasi-odd energies exhibit a monotonic growth and remain mi- 
norized by a positive constant {D/2uj + 1)^/4, 

6[-\d)/A = ik + M + af, k = 0,l,... , (10) 

This split of the spectrum into its three sub-families (H) + (H) + ([10|) is illustrated in 
Table |I] which samples re-orderings of the energy levels in dependence on the growth 
of the coupling D. Once we ignore the superscripted quasi-parities (and, perhaps, a 
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possible separation of the ground state) the spectrum may be visuahsed as an infinite 
sequence of doublets {ek,£i)- Their square roots are in fact equidistant. 

We may summarize that our PT— symmetrized Morse model @ is completely 
characterized by its mapping (^) on the singular harmonic oscillator. One only has 
to note that the mapping need not be invertible easily. Indeed, in the complex 
X— plane, our HO-inspired integration contour C = C(— 1, 1) is by far not unique. 
Firstly, due to the absence of singularities in eq. (|) the non-asymptotic parts of its 
curved integration contour C may be further deformed. One only has to leave the 
asymptotics of C unchanged. Secondly, an asymptotically different curve C{—k, I) of 



ref. [H may be chosen as well |19|. It is characterized by the two odd integers k and I 
and formed by the respective left and right branches x = v — iu with u = ln(c/ cosf) 
and V G (— A;7r/2, — (/c — l)7r/2) \J{{1 — 1)7t/2, /7r/2), pasted together by a (say, straight) 
line with v G {-{k - l)n/2, {I - l)7r/2). 

The feasibility of the non-numerical implementation of the latter idea remains 
an open question. Formally, the new, more general Morse-Schrodinger differential 
equation remains VT invariant and the manifest VT symmetry is also exhibited by 
its boundary conditions. In principle, the new spectrum of energies should remain 
real and bounded below. In practice, the proof must be delivered independently. 
Via a backward mapping, the modified system would be equivalent to a general- 
ized harmonic oscillator with the nontrivial integration curve encircling the essential 
singularity in the origin and moving (perhaps, repeatedly, through the cut) to the 
second (or other) Riemanian sheet (s). 
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Table 1: The D— dependence and pairing of the low-lying spectra (|^). The level- 
crossings occur precisely at the even ratios D / Auj. The energy doublets remain formed 
in their vicinity. The equidistance of the moments \J e\t\D) only takes place at the 
odd integers D/Aoj. 
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